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Abstract 

The non-equilibrium critical-layer analysis of a system of frequency-detuned resonant- 
triads is presented using the generalized scaling of Lee (1997a). It is shown that 
resonant-triads can interact nonlinearly within the common critical layer when their 
(fundamental) Strouhal numbers are different by a factor whose magnitude is of the or- 
der of the growth rate multiplied by the wavenumber of the instability wave. Since the 
growth rates of the instability modes become larger and the critical layers become thicker 
as the instability waves propagate downstream, the frequency-detuned resonant- triads 
that grow independently of each other in the upstream region can interact nonlinearly 
in the later downstream stage. In the final stage of the non-equilibrium critical-layer 
evolution, a wide range of instability waves with the scaled frequencies differing by al- 
most 0(1) can nonlinearly interact. Low-frequency modes are also generated by the 
nonlinear interaction between oblique waves in the critical layer. The system of par- 
tial differential critical-layer equations along with the jump equations are presented 
here. The amplitude equations with their numerical solutions are given in Part 2. The 
nonlinearly generated low-frequency components are also investigated in Part 2. 
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1. Introduction 


Boundary-layer transition usually starts with the spatially growing instability waves as 
reviewed by Kachanov (1994). When the environmental disturbances are relatively small, 
the upstream region is dominated by the linear two-dimensional instability waves. The linear 
two-dimensional behaviour usually persists over long streamwise distances, but the flow 
eventually becomes three-dimensional. The A-shaped structures that have been observed 
in flow- visualization experiments can either be aligned or staggered in alternating rows. 
It is believed that the staggered arrangement is the result of a resonant-triad interaction 
between a pair of subharmonic oblique modes and the fundamental two-dimensional mode 
(Kachanov 1994). As the flow evolves downstream, the initially small side-band instability 
modes become large and the peak in the frequency spectrum becomes wider. Eventually, 
all of the unsteady and steady flow components start to interact nonlinearly and the flow 
becomes turbulent where it will be governed by the fully elliptic type equations rather than 
by the equations that are parabolic in the streamwise direction. 

A self-consistent theoretical analysis of the nonlinear interaction between frequency- 
detuned instability waves during the boundary-layer transition process will be presented 
in this paper. The matched asymptotic method in the high Reynolds number flow will 
be used as in Goldstein & Lee (1992), who studied the resonant-triad interaction in the 
adverse-pressure-gradient boundary layer. Since asymptotic analyses of instability waves in 
high-Reynolds number flows, including critical-layer analyses, are well reviewed by Cowley 
& Wu (1994), we will only summarize the results of the critical-layer analyses that are 
closely related to ours. 
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Goldstein, Durbin & Leib (1987) investigated the nonlinear critical layer in the adverse- 
pressure-gradient boundary layer for an evolving plane wave. They considered a disturbance 
that evolves from a linear instability wave similar to the subsequent free shear layer analyses 
of Goldstein & Leib (1988) and Goldstein & Hultgren (1988). These analyses show that 
the nonlinear effects first become important in the critical layer. The upstream exponential 
growth of the linear instability wave is transformed into algebraic growth in the nonlinear 
region. 

It was found by Goldstein & Choi (1989) that a pair of oblique modes in a free shear 
layer can interact nonlinearly within their common thin critical layer to produce a self- 
interaction term in the amplitude equation. The self-interaction term causes the oblique 
waves to grow very fast and ultimately leads to a singularity at a finite downstream posi- 
tion. This oblique mode interaction can occur at much earlier streamwise position, or at 
smaller amplitudes of the instability waves, than the nonlinear interaction of the plane wave 
considered by Goldstein et al. (1987), Goldstein & Leib (1988) and Goldstein & Hultgren 
(1988). Goldstein & Choi (1989) also showed that the nonlinear interaction between oblique 
modes in the critical layer produces a steady mean flow correction term that is periodic in 
the spanwise direction. The magnitude of this nonlinearly generated component is as big 
as that of the primary oblique modes outside the critical layer. 

In incompressible boundary layer flows, the plane- wave amplitude is usually much larger 
than the oblique-mode amplitude in the upstream region because the plane wave is the most 
unstable mode. Thus, the oblique mode interaction of Goldstein & Choi (1989) is not likely 
to be observed in the early stage of the transition process unless three-dimensional modes 
are artificially excited. The initial nonlinear interaction will then be a parametric-resonance 
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interaction as shown in the critical-layer analysis of the resonant-triad by Goldstein & Lee 
(1992). During the parametric-resonance stage the growth rate of the oblique mode is 
greatly enhanced. The parametric-resonance interaction produces no velocity jump across 
the critical layer at the fundamental frequency. That means the two-dimensional wave 
continues to grow at its initial linear growth rate even when the oblique modes become 
very large relative to the plane wave. Goldstein & Lee (1992) showed that the oblique 
modes eventually react back on the plane wave in the fully-coupled stage. The oblique 
modes also interact nonlinearly between themselves, during this stage, to produce a self- 
interaction term. This term, which was originally obtained by Goldstein & Choi (1989) in 
the free shear layer analysis, causes the instability growth to increase beyond the faster-than- 
exponential growth of the parametric-resonance stage and ultimately leads to a singularity 
at a finite downstream position. Essentially the same amplitude equations were obtained 
by Wu (1992) for the resonant-triad interaction in the oscillating Stokes layer. 

Wu, Lee & Cowley (1993) included the effect of 0(1) viscosity in the critical layer in the 
oblique mode amplitude equation for the oscillating Stokes layer. Increasing the viscosity 
generally delays the occurrence of the finite-time singularity. By taking the highly viscous 
limit of the finite-viscosity amplitude equation, they found that the diffusion layer that 
surrounds the critical layer must be considered in the viscous quasi-equilibrium critical- 
layer analysis. The full finite-viscosity amplitude equations including the back-reaction 
term in the plane- wave amplitude equation were obtained by Wu (1995). 

Mankbadi, Wu & Lee (1993) and Wu (1993) showed that, in the Blasius and favorable- 
pressure-gradient boundary layers, respectively, the dominant nonlinear interactions take 
place in an intermediate diffusion layer that surrounds the critical layer when the viscous 
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effect is large. There is no back-reaction effect on the plane wave which therefore continues 
to grow linearly. 

In the non-equilibrium critical- layer analysis by Goldstein & Lee (1992) and Wu (1992), 
the growth and mean convection (and also viscous) effects can enter the critical-layer dy- 
namics at the same order. However, when the viscous effect becomes large as in Mankbadi 
et al. (1993) and Wu (1993), the mean convection effect balances with the viscous effect 
inside the critical layer, which we refer to as a (viscous) quasi-equilibrium (or viscous- limit) 
critical layer. 

Wundrow, Hultgren & Goldstein (1994) considered the resonant-triad interaction in 
the adverse- pressure-gradient boundary layer when the initial amplitude of the oblique 
mode is small enough so that the plane wave can become strongly nonlinear within its own 
critical layer (as in Goldstein, Durbin & Leib 1987) before it is affected by the oblique 
mode. Goldstein (1994, 1995) extended the analysis of Wundrow et al. (1994) to the 
Blasius boundary layer when the initial oblique-mode amplitude is sufficiently small at the 
start of the parametric resonance. The streamwise evolution of a pair of oblique Tollmien- 
Schlichting waves in the Blasius boundary layer was studied by Wu, Leib & Goldstein 
(1997). 

The Tollmien-Schlichting waves in the boundary layers are initially governed by quasi- 
equilibrium critical-layer dynamics as in Mankbadi et al. (1993), Wu (1993) and Wu et al. 
(1997). However, the results of the analyses by Goldstein (1994, 1995) and Wu et al. (1997) 
and the numerical solutions of Lee (1997a) (who solved the partial differential critical-layer 
equations numerically) indicate that the critical layer in the later downstream stage is 
eventually governed by the non-equilibrium dynamics and the amplitudes of the instability 
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modes are determined by the inviscid or finite- viscosity amplitude equations (Goldstein k 
Lee 1992; Wu 1992, 1995; Wu et al. 1993) or by the system of critical-layer equations (Lee 
1997a) with negligibly small linear growth rates. 

Lee (1997a) developed the generalized scaling which can be applied to both the non- 
equilibrium and quasi-equilibrium critical-layer analyses of zero- and non-zero-pressure- 
gradient boundary layers. He showed that the scalings in the boundary-layer analyses 
mentioned above can be recovered by appropriately choosing the parameters introduced in 
his generalized scaling. The result of the generalized non-equilibrium analysis was presented 
and the system of partial differential critical-layer equations was solved numerically. Lee 
(1997a) showed, by extending the analysis of Wu et al. (1993), that the quasi-equilibrium 
amplitude equations of Mankbadi et al. (1993) and Wu (1993) can also be derived by taking 
the viscous limit of the finite- viscosity amplitude equations. 

The results of the critical-layer analyses mentioned above show that the later down- 
stream stage of the critical-layer evolution is governed by the non-equilibrium dynamics and 
the amplitudes of the instability waves become singular at a finite downstream position. The 
self-interaction between oblique modes is mainly responsible for the explosive growth near 
the singular point. The thickness of the critical layer becomes very large near the singular 
point since it is proportional to the magnitude of the growth rate of the instability mode. 

However, so far, only a single resonant-triad or a pair of oblique modes of the same 
frequency has been considered in the previously mentioned studies. In this paper, a system 
of resonant-triads whose frequencies (of the fundamental waves) are slightly different will be 
analysed. The frequency-detuned resonant-triads can interact between themselves within 
the common critical layer when the magnitude of the frequency difference is the same order 


6 



as that of the growth rates multiplied by the wavenumbers of the instability waves (Lee 
1997b, 1998a). 

The growth rates of the instability modes become larger and the critical layers become 
thicker as the instability waves propagate downstream. The resonant-triads, which initially 
grow independently of each other since their fundamental frequencies are sufficiently de- 
tuned, start to interact nonlinearly in the downstream position where their growth rates 
become large enough to be equal to the magnitude of the phase speed difference (or the 
frequency difference divided by the wavenumber). This implies that most of the unstable 
waves whose scaled frequency differences are nearly 0(1) can interact between themselves 
right before the triple-deck stage where the magnitude of the growth rates and that of the 
wavenumbers are the same. The flow in the triple-deck stage will be governed by the un- 
steady, inviscid, three-dimensional triple-deck equations (Goldstein & Lee 1992, 1993; Wu 
et al. 1997). 

The overall plan of the paper is as follows. The problem is formulated and the gener- 
alized scaling is presented in §2. The non-equilibrium critical-layer analysis is given in §3 
to §7. The solutions in the main boundary layer are given in §3 and those in the viscous 
Stokes layer are presented in §4. In §5 we show that the linear Tollmien solution in the 
inviscid wall layer can be matched onto the solutions of the main boundary layer and of the 
viscous Stokes layer. The derivations of the critical- layer equations and the jump conditions 
are explained in §6. The system of partial differential critical-layer equations and the jump 
equations with appropriate upstream matching conditions are presented in §7. The ampli- 
tude equations which are obtained by solving the system of critical-layer equations in §7 
analytically will be presented in Part 2 (Lee 1998b) together with their numerical solutions. 
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The low-frequency spanwise-periodic components that are generated by the oblique-mode 
interaction in the critical layer will be investigated in Part 2. 

2. Formulation and scaling 

The streamwise, transverse and spanwise coordinates x, y and z are normalized by the 
local boundary-layer thickness A. The mean flow velocity U and the total flow velocity u = 
{it,!?, 1 !*;} are normalized by the local free-stream velocity t/<x>. The pressure p is normalized 
by pUlo and the time t is normalized by A/ {7oo. We suppose that the mean boundary- 
layer flow is nearly two-dimensional and the local Reynolds number Ra(— Uoo A/i/) is 
sufficiently large, so that the unsteady flow is nearly unaffected by boundary-layer growth 
over the region where nonlinear interaction occurs. The base mean flow is nearly parallel 
and develops slowly on the long viscous scale, 

x v = x/Ra- ( 2 - l ) 

The norm aliz ed complex wavenumber a of an unstable mode is small and its imaginary part 
is much smaller than its real part (Drazin & Reid 1981; Smith & Bodonyi 1982; Goldstein 
et al. 1987). 

Following the generalized critical-layer scaling of Lee (1997a), we introduce the wavenum- 
ber parameter a and the growth-rate parameter a r to characterize the small wavenumber 
and the ratio of the small growth rate to the wavenumber, respectively. We are inter- 
ested in the nonlinear interaction of a system of instability waves which is composed of 
2J + 1 resonant-triads (Lee 1997b, 1998a). The jth resonant triad is composed of a sin- 
gle two-dimensional mode of frequency u 3 and wavenumber ay, and a pair of subharmonic 
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oblique modes of frequency u-,/2, streamwise wavenumber nearly equal to oy/2 and spanwise 
wavenumbers ±/?j, where the subscript j can vary from —J to J . The phase speed of the 
plane wave is assumed to be equal to that of the oblique modes of the same resonant-triad 
to the required level of approximation, which is the phase-speed resonance condition for the 
single-resonant-triad interaction. This resonance condition is satisfied when the propaga- 
tion angle of the oblique mode, defined by (5.29), is about 7r/3. It is assumed, in addition, 
for the multi-resonant-triad interaction that the phase speeds (of the plane waves) of the 
resonant-triads are nearly equal to one another. 

The non-equilibrium critical-layer analysis of a system of resonant-triads will be ob- 
tained with the generalized scaling of Lee (1997a). The mean convection effect balances 
with the growth and viscous (in the finite- viscosity case) effects in the non-equilibrium crit- 
ical layer (Goldstein k Lee 1992). The recent studies by Goldstein (1994, 1996), Wu et 
al. (1997) and Lee (1997a) indicate that the later nonlinear critical-layer stage is governed 
by the non-equilibrium dynamics of Goldstein k Lee (1992) even though the earlier stage 
is governed by the quasi-equilibrium dynamics of Mankbadi et al. (1993) and Wu (1993). 
We only present the details of the non-equilibrium critical-layer analysis since the nonlinear 
interactions between frequency-detuned instability waves axe important in the later stages 
of the streamwise evolution of instability modes. However, as will be given in Part 2, we 
can easily obtain the frequency-detuned quasi-equilibrium amplitude equations by taking 
the vis6ous limit (Wu et al. 1993) of the finite-viscosity amplitude equations as shown by 
Lee (1997a) for the single- resonant-triad case. 

The generalized scaling of Lee (1997a) can be written for the non-equilibrium critical- 
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layer analysis as 


a } = a[dt J -\-0{a T % c, = <7[c J + 0(a r )], = cr[c J + 0(<r r )], (}, = ofa, y C] = aY CJ , (2.2) 

xi = a r+1 x, (2.3) 

/i = a T ~ l ft for 1 ^ r ^ 3 and fi = 0(1) for 0 < r ^ 1, (2.4) 

where a characterizes the small wavenumber as mentioned before, Cj and c : are the (nearly 
equal) phase velocities of the two-dimensional and oblique modes respectively, represents 
the spanwise wavenumber of the oblique modes, y C} is the critical level where the base mean- 
flow velocity is equal to the real part of the phase velocity, a } , c } , j3j and Y CJ are order-one 
real constants (which depend on a ) and /z(= R&dP/dx) is the normalized mean pressure 
gradient. The subscript j is used to denote the quantities of the jth resonant triad. The 
real part of c : is only different from that of c 3 by 0 (< 7 r+1 ) as required by the phase-speed 
resonance condition for the single- resonant- triad interaction. When r is equal to 3 and p. 
is an order-one real positive constant, -C 1 characterizes the small adverse pressure 

gradient as in Goldstein & Lee (1992). However, we now allow ft to be non- positive in 
order to extend the scaling to the Blasius and favorable-pressure-gradient boundary layers. 
The magnitude of ft could be smaller than 0(1) in many important cases, i.e. in the 
later downstream stages of the nonlinear evolution of instability waves in boundary layers 
(Wundrow et al. 1994; Goldstein 1994, 1996; Wu et al. 1997; Lee 1997a). 

The scaled Strouhal number s v which is a real constant, of the fundamental two- 
dimensional mode of the jth resonant-triad is given by 

Sj = cr 2 Sj = a 2 QjCj, (2.5) 
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where the frequency or unsealed Strouhal number s 2 of the plane wave is equal to oc^Cj (see 
(3.7)). We are interested in the case when the scaled Strouhal number of the jth resonant- 
triad is different from that of the ‘reference’ Oth resonant-triad by a small factor whose 
magnitude is of the order of the growth-rate parameter <x r , 

Sj = s 0 (l + o t jx) for (2.6) 

where the frequency-detuning factor x is an order-one real constant. It will be shown 
by (5.22) that the scaled phase velocities c 2 for j = — J, ...,«/ axe different one another by 
0(<r r ), therefore, the additional resonance condition for the multi-resonant-triad interaction 
is satisfied. 

If the Strouhal number s 2 and the small parameter a are chosen (these parameters 
can be chosen from experiments and appropriate scalings for the specific problem), the 
wavenumber a 3 can be determined from the long- wavelength small-growth-rate dispersion 
relations (5.8) - (5.13) in §5. The phase speed c 2 can be obtained from (2.5) and then the 
critical level y CJ which is defined as the transverse location where 

u(y = y cj) = at,, (2.7) 

is determined if the base mean flow U is known. 

The non-equilibrium critical-layer analysis by Goldstein & Lee (1992) shows that the 
nonlinear interaction between instability waves of a resonant-triad first occurs within the 
critical layer whose thickness is of the order of cr r+1 . When s 2 is given by (2.6) the critical 
level of the jth resonant-triad, y CJ is only different from that of the Oth resonant-triad, yco by 
0(a T+1 ) (see (5.22)). Therefore, all of the instability waves of the detuned resonant-triads 
will then nonlinearly interact between themselves within the overlapped critical layers. 
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The viscous effects will enter the critical-layer momentum equation while making only 
relatively insignificant modifications to the external flow when the viscous parameter (Ben- 
ney k Bergeron 1969; Haberman 1972) 

A = l/(<r 3r+4 i? A ) (2.8) 

is order one. The scalings for the amplitudes are, following the previous studies by Goldstein 
& Lee (1992, 1993) and Lee (1997a), for all j, 

e 2 d = <t 4 ' +1 , 6 3d = 602 = a 3r+1 , (2.9) 

where e 2 <i and fad are the amplitude scalings of the two-dimensional and oblique waves. The 
nonlinear interaction between oblique modes in the critical layer produces a term which is 
periodic in the spanwise direction (Goldstein k Choi 1989). This term is large in the sense 
that its magnitude fa 3 is equal to 6 3d . This nonlinearly generated term becomes a function 
of a slow time, as will be shown in (5.30), in the multi-frequency and multi-resonant-triad 
cases considered in this study, though it is independent of time in the single-frequency 
(Goldstein k Choi 1989) or single- resonant-triad (Goldstein k Lee 1992) cases. 

The nonlinear interactions between instability waves occur within the critical layer 
whose thickness is of the same order of magnitude as the growth rate, i.e. 0(a r+1 ). We 
do not need to include the diffusion layer (c.f. Mankbadi et al. 1993; Wu 1993) since it is 
inseparable from the critical layer in the non-equilibrium critical-layer scaling (Lee 1997a). 
The viscous Stokes layer does not make a leading-order contribution to the growth rate in 
the non-equilibrium analysis. However, in order to make the analysis valid in the viscous 
li mi t, we will include the effect of the viscous Stokes layer whose thickness is of the order 
of <r 3r / 2+1 . A schematic diagram of the multi-layer structure is given in figure 1 with the 
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critical levels of the Oth, 7th and — 7th resonant-triads indicated. The distance from the 
wall to the jth critical level is y CJ that is of 0(a) (see (5.22) for their relative positions). 

Finally, we will assume that 0 < r ^ 3 for the simplicity. This range of r covers most of 
the important types of boundary-layer flows as was summarized in table 1 of Lee (1997a). 

3. Main boundary layer 

Since the unsteady flow outside the critical layer is mostly governed by linear dynamics 
to the required order of approximation, the velocity field of a system of 27+ 1 resonant-triads 

in the main boundary layer, where y is of order one, can be written as 

J J 

u = U(y,x v ) + e 2 d X + S 3d X R e 2A,C/j(y, x x )e‘* j/2 cos Z } 

J=-J J=~J 

+^02X^Retfo,2y,r(y,*i,M)+ ...» (3.1) 

: e 

j 

v = V(y,x v ) - e 2 d X Rei ajBj(ii)$j(y,a;i)e ,Xj 

j=-J 

J 

-63 d X R e2i-fjA J (xi)^ J (y,Xi)e lX3/2 cosZ ] + (3.2) 
i=-J 
J 

w = 63d X R e2i A» X\)e lX} l 2 sin Z 3 + . . . , (3.3) 

J=-J 

J J 

p=P + e 2d X ^B ] P 3 (y,x 1 )e i ^ + 6 3d X Re2A J P,(y,z 1 )e i *>/ 2 cosZ J + ..., (3.4) 

]=-J j=-J 

with 

Xj = aotjX — a^Sjt, Z 3 = afijZ, (3-5) 

where Re denotes the real part, x v and ii are defined in (2.1) and (2.3), and e 2d , fad and 
6 02 are given in (2.9). Since s 3 , a 3 and /3j are real X 3 and Z 3 are real, but we allow the 
oblique-mode amplitude A } and the plane- wave amplitude B : of the jth resonant-triad to be 
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complex. The mean flow pressure P is a function of x in general and the pressure gradient 
dP/dx that is equal to n/R& is scaled by (2.4). The last unsteady term in (3.1) is induced 
by nonlinear effects in the critical layer. We are interested in the case when the amplitudes 
of the oblique modes with positive and negative propagation angles are the same. The pair 
of oblique modes with the same amplitude forms a standing three-dimensional wave that 
propagates only in the stream wise direction. 

The mode shapes and $ 3 satisfy the Rayleigh’s equations to the required levels of 
approximation 

(U - c 3 )(D 2 - 7 j)$_, - U'% = 0, (3.6) 


and the complex wavenumbers a 3 and 7 ^ and phase speeds c 3 and c 3 are given to the required 
order of approximation by 


7; = a \% + a 


,r ML' 


ere. 


c, = 


l + a r kA' 3 /(ia 3 A 3 y 


(3.7) 


where the prime denotes differentiation with respect to the relevant argument, D = d/dy, 
and {A 3 , $ 3 , c 3 , 7 j , 7,, k} can denote either {B 3 , $ 3 , c 3 , a 3 , a 3 , 1} for the two-dimensional wave 
or {A 3 , ^j,Cj, 7 j, 7 j, 2 } for the oblique modes. We define 7 ^ as 


[K/2) ! + ^] ,/! . 


(3.8) 


Instead of solving (3.6) with the leading-order inviscid wall boundary condition (Goldstein 
& Lee 1992; Wundrow et al. 1994; Lee 1997a), we will appropriately match the solution of 
(3.6) with that of the viscous Stokes layer. 

By substituting (3.1) - (3.4) into the continuity and momentum equations we can show 
that (Goldstein & Lee 1992) U 3 , W 3 , P 3 and P 3 can be expressed in terms of D$ 3 , $ 3 
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(3-9) 


and as 

3 " [ J « 2 (tf ~ c,)J 27, + [ J V - c, ^ + &) ) J 7 j^ 

W, = (U 3 - n iD^j/aj) otj/fij - <r r (&, - \a 3 D* 3 /%) A'JifoAj), (3.10) 

P, = 7^%/(7, 2 - #) 1/2 - “ c 3 ) U 3 , (3.11) 

P 3 = U'9j-(U -c 3 )D9 3 . (3.12) 

When r > 1, the mean flow U ( y , x v ) can be written as (Goldstein & Lee 1992; Mankbadi 
et al. 1993; Wundrow et al. 1994; Lee 1997a) 


U = U B + v r - 1 U P + U(y,x v ) for y = 0(1), 


(3.13) 


where U B is the Blasius velocity and a r ~ 1 U p is a small correction due to the pressure 
gradient. As y approaches the wall, which is assumed to be located at y = 0, 


U B 


+ + e *-' i (- f i '+? 2 + •••)■ 


(3.14) 


where the constant r 0 (~ 0.332) denotes the scaled Blasius skin friction and p, was defined 
in (2.4). The total wall-shear stress r w that is the sum of the Blasius skin friction and the 
correction due to the mean pressure gradient can be written as (Wundrow et al. 1994) 


T W = T 0 - a T 1 pT. 


(3.15) 


When 0 < r ^ 1, the correction term due to the pressure gradient is not small any 
more since the magnitude of the mean pressure gradient /x could be of order one as in (2.4). 
Therefore, (3.14) and (3.15) should be replaced by (Wu 1993) 

U -*■ r w y + 5 /xj/ 2 + . . . as y -* 0. (3.16) 
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In the following analysis p will be treated as an order-one quantity. However, the results 
will also be valid for 0 < r ^ 1 if p is replaced by <r 1 ~ r ^i. 

The terms U(y,x v ) and V(y,x v ) in (3.13) and (3.2), which are functions of the long 
viscous scale x v given by (2.1), are included due to the slow viscous spreading of the base 
mean flow and they can be found by expanding the mean flow solution about x = 0. 
However, they do not make any major contribution to the unsteady flow over the entire 
streamwise region in which the nonlinear interaction takes place since the local Reynolds 
number is assumed to be sufficiently large. 

As in Goldstein et al. (1987), Goldstein k Lee (1992), Wundrow et al. (1994), Lee 
(1997a), and others, the matching can be simplified by using the Miles’ (1962) solution, 
which is obtained by transforming (3.6) into a Riccati equation, and the classical inviscid 
function (Lin 1955). The latter is written as 






(3.17) 


U'$j — (U — Cj)D$j 

where VVj can denote either VV, for the plane wave or W 3 for the oblique modes. Substituting 
the solution by Miles (1962) along with (3.7) into (3.17), expanding for small <r, and using 
(3.13) to (3.16) as in Goldstein k Lee (1992), we obtain 

>Vj = - Cjf/ _ T lj x + a 2c : J 2 + a 2 c 2 J 3 -f + ct 2 %c 3 t w (2 J 4 + ac 2 J h ) 

7,(1 -°c 3 Y (1 — aCj) A \ 8 t w J 

'h (3 - 18) 


—a 


where the constants J\ to Js are given in Appendix A of Lee (1997a) and we have put 


3— r 

Pcj = P ~ ° 



(3.19) 
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We have also used the fact that 


c, = r w Y„ + a’i Y e ] (a - S-'tX/*') +■■■• (3.20) 

which is obtained from (2.2), (2.7) and (3.13) - (3.16). The above equation (3.20) determines 
Y CJ once c } is known. 

4. The viscous Stokes layer 

As shown in Goldstein & Lee (1992), Wundrow et al., (1994) and Lee (1997a), the 
viscous Stokes layer adjacent to the wall does not make a leading-order contribution to 
the growth of the instability wave in the non-equilibrium critical-layer analysis. This thin 
viscosity dominated layer, however, produces an order-one growth term when the Tollmien- 
Schlichting waves are considered as in the quasi-equilibrium analyses by Mankbadi et al. 
(1993) and Wu (1993). It was pointed out by Wu et al. (1993) and shown in detail by 
Lee (1997a) that the quasi-equilibrium amplitude equations of Mankbadi et al. (1993) and 
Wu (1993) can be alternatively derived by taking the viscous limit of the finite- viscosity 
amplitude equations provided the viscous Stokes layer effect is appropriately included in the 
the linear growth rates. Goldstein & Lee (1992, 1993) identified these correction terms by 
comparing their results with those of the quasi-equilibrium critical-layer analyses. In this 
section, the unsteady flow in the viscous Stokes layer will be analyzed in order to obtain, 
in a systematic way, the linear growth rate of the instability wave that is also valid in the 
viscous limit. 

The generalized scaling by Lee (1997a) shows that the viscous Stokes layer for the 
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non-equilibrium critical-layer analysis is scaled as 


y = a 3r/2+1 y. (4.1) 

If we substitute (3.1) - (3.4), (4.1) and 

= a 3r/2+1 l„ (4.2) 

into the continuity and momentum equations, we can obtain 

+ x (*, - j) = o, (4.3) 

where {$ J7 $j,k} can denote either 1} for the plane wave or 2} for the 

oblique modes and the base mean flow U is given by (3.13) to (3.16). The bounded solution 
which satisfies the no-slip boundary condition at the wall ( y = 0) is 

4, = ( 7 -(3’72+Dtf + (rjbj) (e~ b >* - l) , (4.4) 

where 

1/^ = (1 + i) (ifcA/^) 1/2 . (4.5) 

It is easy to show from (4.4) that 

— *• £T - ( 3r / 2+1 ){7 — r w /bj as y — ► oo. (4.6) 

This will be matched with the inviscid wall layer solution in the following section. 

5. The inviscid wall layer 

The analyses of Graebel (1966) and Nield (1972) suggest that we have to introduce the 
scaled transverse coordinate 

Y = y/<T (5.1) 
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directly into (3.6) before attempting to obtain the solution. The solution is of the form 

9 } = <tt u Y + o t+ 1 F 3 (Y, ]>f ) + ..., (5.2) 

where the integration constant 4>f, which can denote either <t>f or <f>f for the plane or oblique 
waves respectively, is a complex function of x\ in general. Matching (5.2) with the viscous 
Stokes-layer solution (4.6), along with (3.13) - (3.16), (4.1), (4.2) and (5.1), shows that F 3 
must satisfy 

F 3 -e r ' 2 T w /b 3 as F — ► 0, (5.3) 

where b 3 is given by (4.5). 

We can show that F which satisfies the above matching condition and can be discon- 
tinuous across the critical level Y cv is given by 

Ff = f 3 (Y) + \fi c: Y C] [( Y - Y CJ )$f + Y e] 4>j] - a^rjb, + . . ., for F>F CJ , (5.4) 

where 

f 3 (Y) = fi CJ [i Y 2 + Y CJ {(Y - Y ej ) In \Y - Y CJ | + Y C1 In F CJ }] - a 3 ~ r r 2 0 Y 3 ( Y CJ + \y) /4!, 

(5.5) 

and ji C3 is given by (3.19). The last term in (5.4) was added due to the effect of the viscous 
Stokes layer and becomes (1 + i)r w (\kc r \/s^ if we substitute (4.5). Its magnitude is 
of higher order, i.e. 0(a r / 2 ), in the present finite-viscosity case, but it will be of order one 
in the viscous-limit (quasi-equilibrium critical-layer) case where A becomes 0(cr~ r ). The 
rescaling of A in the viscous limit will be given in Part 2 (also by (8.8) of Lee (1997a)). It is 
obvious that we need to include more (inviscid) terms (Wundrow et al. 1994) on the right 
hand side of (5.4) in order to make that to be uniformly valid up to 0(cr T / 2 ). However, in 
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the following analysis, we will keep the viscous correction term as in (5.4) and treat it as 
if the magnitude of a T A is of order one since we are mostly interested in the leading-order 
effect of that term in the viscous limit. 

In the later non-equilibrium critical-layer stage, r becomes smaller than 3 as summa- 
rized in table 1 of Lee (1997a) and jl is also likely to become negligibly small a s in Wundrow 
et al (1994). Thus, /2 CJ defined in (3.19) becomes negligibly small so does the linear velocity 
jump across the critical layer. When /i CJ < 1, (5-4) can be explicitly rewritten as 

Ff = fj(Y) + i/2 CJ Y CJ \(Y - Y CJ ) ($">* + »cA L)± ) + Ycj $ N) - + M L) ')] 

-a T, 2 T w fb 3 + • . • , for Y <Y CJ , (5.6) 

where accounts for the 0(1) nonlinear jump while p. C j4>j L ^ ± accounts for the very 

small linear jump. However, for the simplicity, we will treat p, CJ as an 0(1) quantity and 
(5.4) will be used instead of (5.6) in the following analysis. The result will be still valid 
even when fl CJ < 1 if p , is simply replaced by 4>\ N ^ +/Z CJ $^ ± . 

Inserting (3.13) - (3.16), (5.1), (5.2), (5.4) and (5.5) into (3.17) and re-expanding, we 


obtain 



(5.7) 

for Y > Y C} . The inviscid function (5.7) in the inviscid wall layer will be matched with that 
in (3.18) obtained in the main boundary layer as explained below. 

The (linear) real part of the wavenumber of the jth resonant-triad a 3 in (2.2) possesses 
a power- series expansion 

dj = dj°) + aa.^ + a 2 a^ + . . . + cr r (ln cx)a^ rL ^ + cr r a^ + (5.8) 
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Matching (5.7) with (3.18) along with (2.5), (3.13) - (3.16), (4.5) and (5.1) for the jth plane 
wave where {Aj,*t r 4>f,k} must be replaced by {Bj,a v <i>f, 1} shows that (Goldstein & Lee 


1992; Wundrow et al. 1994; Lee 1997a), for 0 < r ^ 3, 

a<°> = (r (5.9) 
oj 1 * = s,(l - if r > 1, (5.10) 

*S 2) = -i^(^) 1/2 [3Ji + 2J 2 - \r w {Jl + 8J4)] if r>2, (5.11) 

"j rL) = 2 Ph/ T l - I (5.12) 


«('» = 1(1 + In Y',)PS,/tI - i (a,/c,) 2 (\S A/S ,) - o»-*\s) [(3 + 2r 2 J t ) h + 

+* - t w J 5 + {7/48 + } (t,Y c ,[c ,) 2 In Y cl } /r *] + cr^a^ + at’-'lfiW, (5.13) 


and 




(5.14) 


where <rt n— is defined to be equal to <r n-r if n ^ r or equal to zero when n < r, fi CJ is given 


by (3.19) and J\ to J$ are given in Appendix A of Lee (1997a). 

The above equations (5.9) - (5.13) are the usual long- wavelength small-growth-rate 
dispersion relations that determine the wavenumber a, in terms of s r The phase speed 
Cj which possesses a similar expansion to (5.8) is determined by (2.5). The complex plane 
wave amplitude B : of the j resonant-triad is determined by the phase-jump equation (5.14). 
The effect of the viscosity now appears as the second term on the right hand side of (5.14) 
and always increases the growth rate. The phase jump <j>J — <f>f across the critical layer 
can be obtained by considering the critical layer and is equal to the sum of —it and a 
complex function of Xi as will be shown in §6 (Drazin & Reid 1981; Goldstein & Lee 1992). 
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The former and the real part of the latter represent the linear and nonlinear growth rates, 
respectively, while the nonlinear correction of the wavenumber is given by the imaginary 
part of the latter. The linear part of the wavenumber, a } has been fully determined, up to 
the order of a T , by matching as given by (5.8) - (5.13) (Wundrow et al. 1994). 

As shown above, a, and c, have been determined from the matching of the plane-wave 
inviscid functions. The wavenumber 7 , can now be determined by matching the inviscid 
functions of the oblique modes. It is convenient to put 

/ 3 J — ^(a, + (T r K s ), (5.15) 


and then from (3.8) 


“1“ G 4^3 "4* • ♦ • * 


(5.16) 


where k s will be determined from matching. By matching (5.7) with (3.18) along with (2.5), 
(3.13)- (3.16), ( 4 . 5 ), ( 5 . 1 ) and (5.16) and replacing by {A^ for the 


oblique modes, we can show that 


(2l + ?2i\ = J-l \h&U- -tf) - 
+ s i J K 4 T - \ s > 




(5.17) 


where k vj appears due to the viscous Stokes-layer effect. Since we have chosen k 3 in (5.15) 


to be the same for all j, which we put 


k , = > (75 - 0 KV *.) 1 ' 2 , ( 5 - 18 ) 

the higher-order j-dependent correction term k vj must be included in (5.17), 

k vj = - l) (* r A ) 1/2 (%/«,) [a, 7 ,/(c 2 ^ /2 ) - ag/(c& /2 )] . (5.19) 

Using ( 2 . 6 ), (5.16), ( 5 . 21 ) and ( 5 . 22 ), we can show that the magnitude of is of the order 
of <r r (a T A) 1/2 . The phase jump <f>J - <j>+ is equal to -x (for the linear growth rate) plus 
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a complex function of x x (whose real and imaginary parts represent the nonlinear growth 
rate and the nonlinear wavenumber correction, respectively). 


It is easy to show from (3.1) and (3.5) that the nonlinearly corrected streamwise 
wavenumber of the oblique mode becomes <7dy/2 — o^^ReiAj / Ay Thus, the nonlinear 
propagation angle of the ;th oblique modes, ±0 3 can be computed, up to 0(a r ), by 

L = tan' 1 ( Si ^-J) , (5.20) 

J Vdj/2 - a r ReiA' 3 /Aj 

where and p 3 can be evaluated from (5.8) - (5.13), (5.15) and (5.18). As we expected 
(Goldstein & Lee 1992), 0 3 becomes x/3 at the leading order. 

The result presented so far is basically a linear superposition of the (modified) results 
of the single- resonant-triad analyses by Goldstein & Lee (1992) and Lee (1997a) since the 
unsteady flows in the main boundary layer, viscous Stokes layer and inviscid wall layer are 
governed by the linear dynamics to the required order of approximation. 

From (2.5), (3.19), (3.20), (5.8) - (5.13), (5.15), (5.16) and (5.18), the quantities of the 
jth resonant-triad can be expressed in terms of the ‘reference’ Oth resonant-triad, when s } 
is given by (2.6), 


a, = a(l + <r r ljx), fa = p + a^jxa, 7j = 7 + ° T \3X& + • • (5-21) 

c J /c=l + <7 r ij X + ..., Y C] /Y c = 1 + <t t Ij X + ..-, fie] = fic + 0(cr 3 ), (5.22) 

where 

x = X (l + ^ (1) /a + <T 2 2d^/a + ...), (5.23) 

and aW and are given by (5.10) and (5.11). We have omitted in the above equations 
and we will omit in the rest of this paper the subscript 0 for the quantities of the ‘reference’ 
Oth resonant- triad for notational simplicity. 
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It is convenient to introduce the slow time and the slow spanwise coordinate z\ which 
are defined as 

*i = a r+2 xst, Z\ = CT r+1 ^x«z. (5.24) 

By using (2.6), (5.21) and (5.24), we can show that X 3 and Z 3 defined in (3.5) become 

Xj = X - jh + \jxax i , Z 3 = Z + jzi . (5.25) 


If we put 


A,(n) = A 3 ( Xl y* s *''\ B } ( Xl ) = b 3 ( Xi )^ s ^ 2 , 


(5.26) 


and use (5.19), (5.21) and (5.22), the phase-jump equations (5.17) and (5.14) can be rewrit- 
ten as, for —J ^ j ^ J, 


( 


cos 6 + 


1 


'A!, i 


-j i 7 2 

j'H * -w 


¥(♦;- «)-(¥ 


a r Ay/ 2 ' 


a 2 


B' } i .. 

~ 2 JXa = " 2 ^ 


«-©1' 

where p. c for the Oth resonant-triad is given by (3.19) and 6 is defined as 


(5.27) 

(5.28) 


6 = cos 1 (50/7) , (5.29) 

which is the same, but cr r Reii4'/ A 3 , as the nonlinear propagation angle 9 in (5.20) when 
j = 0. The relation (5.16) shows that 9 is 7r/3 at the leading order, but it can be fully 
determined because a and 7 are known from (5.8) to (5.13), (5.16) and (5.18). 

The phase jumps across the critical layer will be determined from the nonlinear analysis 
within the critical layer. Since we are interested in the instability waves that are spatially 
growing in the streamwise direction in the linear upstream region, they will continue to 
grow only in the streamwise direction as they propagate downstream. The amplitudes A 3 
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and Bj are complex functions of x\ but are independent of t\ and Z\. The wavetrains (Wu, 
Stewart & Cowley 1996) which can grow in the slow time and in the spanwise direction will 
be considered in a subsequent paper. 

Equations (5.1), (5.2), (5.4), (5.5), (5.21) - (5.23) along with (3.9) - (3.16) show that 


(3.1) - (3.4) can be rewritten in this inviscid Tollmien region as (Goldstein & Lee 1992; Lee 
1997a), 


2 J 

u = ar w Y + a r+1 i/2Y 2 - a 4 ^y Y 4 + e 2 d £ Re [r w + <r r (/ y + i/2 c Y c 4>f)\ B ] E : 
+^3d £ 2r tu (sec0)Re jl + C(sin 2 0) [l + 0- r |jx(l + C)] + o" r ^(sin 2 6) j / + — f^n 2 q 

o - - s (- * 12) * ¥ (*: * <&) 

— (1 4- i)r w ( —~j \A 3 E cj + 6 02 £ £ Re^o. 2 ;^-^ 2 cos[2Z + (j + 2l)z x ) 

' S ' J J j=-2 J ( 

+ ..., (5.30) 


v = -a 2 ^2 Re \t 2 d {r„,Y + <r r ificYc (Y - Y c ) 4>f + . . .} iaBjEj 
i=-J 

+6 3d {t w Y + a T ip. c Y c (Y - Y c ) d>f + . . .} 2ry A 3 E CJ ] + . . . , (5.31) 

j 

w = -6 3 d £ 2(sin0)Reir u ,O4 J £ aj + ..., (5.32) 


}=-J 


p = P + £ ffcr^Re 


tuBjE) + 6 3( f2(cos 9)AjE c ^ + . . . , 


(5.33) 


j=-J 

where e 2 d, hd and <5 02 are given by (2.9), / and /2 C are given by (5.5) and (3.19) and we 
have put 


Ej = # x -*\ E CJ = e* x ~*'V 2 cos (Z + jz 1 ), E SJ = e^ x ~^ 2 sm(Z + jz x ), (5.34) 


and 


C = c/(r w Y - c). 


(5.35) 
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In the streamwise velocity (5.30), the term a T ^jx(l + Q appears because the quantities 
of the jth resonant-triad are expanded about those of the 0th resonant-triad as in (5.21) - 
(5.23) and the term which is proportional to (cr r A/s) 1 / 2 is included due to the viscous Stokes 
layer effect. The last term on the right hand side of (5.30) that is a function of the slow 
time t\ is included in order to match with the critical-layer solution. The 0(<7 r+2 ) terms 
that are continuous across the critical layer are not shown in (5.31) since they do not play 
any major role in the analysis. The mean pressure is denoted by P in (5.33). The above 
solutions have to be rescaled in the critical layer since they become singular at Y = Y c . 

Although we have used the non-equilibrium critical-layer scaling given in §2, the results 
in this section as well as in §3 and §4 are also valid in the viscous limit since the leading-order 
viscosity effect was included by considering the viscous Stokes layer. In the viscous limit 
where A defined in (2.8) is of 0(a~ T ), see (8.8) of Lee (1997a), the critical layer becomes of 
the quasi-equilibrium type and the viscous Stokes layer becomes thicker (Mankbadi et al. 
1993; Wu 1993). The generalized scaling given in table 1 and (8.9) of Lee (1997a) indicates 
that A, €2 d, $ 3 d, ft 02 and y in (2.8), (2.9) and (4.1) must be rescaled for the quasi-equilibrium 
critical-layer analysis. 

6. The critical layer 

The nonlinear analysis of the detuned resonant- triads within the critical layer will be 
presented in this section. The system of critical-layer equations and jump equations will be 
derived here and their normalized version will be given in the following section. 

The thickness of the non-equilibrium critical layer is of the same order as the growth 
rate, i.e. 0(a r+1 ), as shown by Goldstein et al. (1987), Goldstein & Lee (1992), Wundrow 
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et al. (1994) and Lee (1997a). The growth rate and the linear convection terms become of 
the same order of magnitude there. The appropriate transverse coordinate in this region is 

f) = (Y -Y c )/a T = (y-y c )/a T+1 . (6.1) 


The continuity and momentum equations expressed in terms of the scaled variables x\, 
X, Z, z u h and fj , which were defined in (2.3), (3.5), (5.24) and (6.1), become 


au x + v n + fiw z + cr r (u x j + ^-x aw zi) = 

Z){u,t>,to} = - + ° T Vx x ,a _ * 2r+4 *p < p/?Pz + ° T ^X«P*i} > 


where we have put 
D = a(u - a c) 


w *w +' Sw w + ' '' ("air + + ■ 7,+1 { x& ~ c k ~ A 


( 6 . 2 ) 


(6.3) 


&)' 

(6.4) 


v = <r r+2 i ;. 


(6.5) 


The viscous parameter A is defined by (2.8), the frequency-detuning factor x is introduced 
by (2.6), x is given by (5.23) and u,v,w and p are total velocities and pressure. The 
quantities a, /?, c and Y c represent those of the Oth resonant-triad, i.e. when j = 0, as 
mentioned in the previous section. 

Putting (2.9) and (6.1) into (5.30) to (5.33) and re-expanding the results show that the 
critical-layer solution must be of the form 

u-ac = a T+X T w fj + a 2r+l u^ + <r 3r+1 u< 2 ) + <r 4r+ V 3 > + <r 5r+ V 4 > + . . . , (6.6) 

v = —<r 2r+1 £ 2jr w Y c ReiAjE CJ + <r 3r+ V 2) + <r 4r+1 u (3) + a 5r+1 v (4) + . . . , (6.7) 

i=-J 

w = a 2r+1 w^ + a 3r+ V 2 > + a 4r+ V 3 > + <r 5r+ V 4 ) + . . . , (6.8) 
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J 

p = P + <r 3r+2 53 2cr w (cos 0)Re i,.E ej + a 4r+ V 2) + * 5r+ V 3) + . . . , (6.9) 

j=-J 

where E C] is defined by (5.34), t w is given by (3.15) and 6 is defined in (5.29). The solutions 
in the inviscid Tollmien region (5.30) - (5.33) indicate that we must include the terms 
t7 4r+ i (in a j u (3L) an( j < 7 S r +i(i n a)u( 4i ) into (6.6) and similar terms for v, w and p into (6.7) - 
(6.9). However, these logarithmic terms do not play any active role in the present analysis 
and thus they are omitted for the simplicity. Substituting these expansions (6.6) - (6.9) 
into (6.2) - (6.5) with (3.15) we can obtain the governing equations for u^ l \ v^ l \ wW and 
pW which are given in Appendix A. 

As in previous studies (Goldstein & Lee 1992; Lee 1997a) we can solve these equations 
(A 1) - (A 11) subject to the transverse boundary conditions that they match onto the outer 
solutions (5.30) - (5.33). The relevant solutions for the lowest-order are given by 

J 

u (1) = u$(t},xi)+ 53 Re2i(tan Q)Q^\fj, x\)E CJ , (6.10) 

]=-J 

J 

io (1) = 53 Re2Qj 1) (^,xi)£: aj , (6.11) 

j=-J 

where E SJ and E CJ are given by (5.34). From (5.30) and (6.6), we can show that the base 
mean flow component becomes 

U M = ( 6 - 12 ) 

where 

a M = Wc ~ ^ 3 - r r 2 y c 3 /12 + . . . . (6.13) 

However, since the constant a $ is not required to be evaluated in the following analysis, 
we will retain the general form (6.12) without substituting (6.13). 
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The nonlinear equations in Appendix A indicate that more Fourier components must 
be considered in this multi-frequency analysis than in the single-resonant-triad case. We 


can show that u (2 \ u (3) and possess the following expansions: 

« (2) = U M + E Re { h,])' El,] + Re { (^S/O;,) Eq-j + 2y 

J|1 1 2 J|2 

E 2 -j + Y2 E 2 -j + ^ (Uo,i,j,(Eo-j + U-0, 2 j,(Eo -j) ^2;j+2<| , (6.14) 


« (3) = + E ^ + £ £ Re { *»+* 

+ ) 1 v} + £ Rei ^ ' + E 

[ m:-l m:l X 

^l;-j+2<-2m) t )} + £ ( (^Ua/ 1 ^) + E (^lfl;i/,m^l;-J+2m) j + £ £ Re 

/:[— 2] l m:-l J J j|3 «:[2] 

£ 3y (^3ay,/^l;-J+2/) + E { (^34uAm^ly-2m) + (^3,i;j,; im ^l;j-2<+2m) } 1 + • • • > (6.15) 

. m:l ^ -J 


««> = . 5 ? + E + EE Re + E R * E { «j/ 

j|l j|2 fcl J|3 Up] 

^’twl* E'i 'ti-t + l ^2;j+m)+ E (^2,0;h^<v) ^2;2 j+2 + E Re ^2;j 


E E ((^ilta^-2™)' + (^, »«*-»«-)’) + E 

[*2 [m:l ' ' m; — 1 J 

+ E E +EE^[(«A-»)'+ E {«U. 

l\—2 m:— 1 J M * ; [3] L m P) 

£0 i;-2<+2ra) t + ^ (^ 2 ^i<,m,n^ 0 ;j- 2 m- 2 n) + E ( (^ 2 , 0 ;j,r,m,n^Oy- 2 n) + (^2,0y/,m,n 

n:-l n:l V 


£o;j— 2 m+ 2 n 


)*)}+ E { (^2,0;],«,m^0y-2<-2m) + E (^2,0;j/,m,n^Oj-2n) | +•••» (6.16) 

m:[- 2] ( n: “ 1 J . 
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where we have put 


E a , } = e ^ aX ~^ 2 , £ b -, t = e i{bZ+izi) , (6.17) 

(fj (k ) t . ,£ C S = U {k ) r, ^£c; + U {k \ i, /-c-, for k = 2L,...,4o. (6.18) 

It is easy to show that the previously defined E 3 , E CJ and Es: by (5.34) are equal to 
\Ei ;j (£i ;: + £-i ; -j) and -\\E x - : (£\- 3 - £-\--j), respectively. The summation notations used 
in the above equations are defined as, for a = 1,2,3,..., 


E s E • («•«> 

j\d j= — dJ 

E = E . E = E ■ («- 2 °) 

{/,m,n}:±d {<,m,n}=max(— aJ,— aJ±{j,<,m}) {f,m}:[±d] {/,m}=max(— dJ,- J±{j/}) 

The streamwise velocity expansions in (6.14) - (6.16) are basically composed of the 
Fourier components m n ]E a -,i£b;l where the integers a and b are coefficients of ^iX 

and i Z as in (6.17) and i and l are expressed in terms of j, £, m and n. All daggered terms 
are composed of two components, for example, (E^joy, -® 2 y in (6-14) is composed of 
the two components U 2 %j,(E 2 ;j£o;j- 2 t and U^_ 0 .. ! t E 2 ^£Q-- 1 + 2 i- The shape function U^l 0 y/ 
is different from ( in general, although they will be shown to be the same in the present 
case by (6.26). By putting J = 0, or consequently j = l = m = n = 0, we can recover the 
expansions for the single resonant-triad case given by Goldstein & Lee (1992), Wu (1992) 
and Lee (1997a). The oblique- mode velocity jump across the critical layer can be obtained 
from those of the nine Fourier components given in (6.15) with the subscript 1, 1; 
and the superscripts 3 L, 3a, . . ., 3/t. Meanwhile, we need to consider all sixteen components 
in (6.16), denoted with the superscripts 4L, 4a, . . ., 4o, in order to obtain the velocity jump 
of the plane wave. 
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«, ( 2 ) 

The solution of the critical-layer equation for the UQ 2jt (Eo\j£2;)+2t com P onent i given 
in Part 2, shows that the nonlinear interaction within the critical layer produces a non-zero 
jump of this component (Goldstein & Choi 1989). It is worth to note that this term is a 
function of the slow time ti defined by (5.24). For the single-resonant-triad interaction, the 
subscripts j and t become zero and thus this component becomes steady as we expected 
from the results of Goldstein & Choi (1989), Goldstein & Lee (1992, 1993) and Wu (1992). 

By substituting (6.1) into the base mean-flow component in (5.30), or directly into 

( 2 ) 

(3.13) to (3.16), re-expanding the result and comparing with (6.6), we can show that u M 
and uff in (6.14) and (6.15) are given by 

U M = a lA/*? 2 + a 2MV + a 3M x l + • • • > ( 6 - 21 ) 

+ 4 3 A^i + ■ ■ • » ( 6 - 22 ) 


where 


a \M = §& + ••■. 


(6.23) 


and ji c is given by (3.19). The constant afy is nearly equal to one half of the double 
derivative of the base mean flow at the critical level y — oY c divided by the scaling factor 
<r r_1 . The other constants a^, a^, a^, and can also be found although the 
results were not presented since they do not play any major role in the analysis. The last 
terms in (6.21) and (6.22) are functions of x\ and they are included in order to match the 
slight viscous correction that enters the external solution due to the slow viscous spreading 
of the mean flow as in (3.13). 

The expansions for the other velocities and and the pressure pW for k =2, 3 
and 4 can be obtained if we replace in (6.14) to (6.16) by V^l j[Am n] , 
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and Patj[ ( m n] and replace the mean flow component by vff, w$ and p$, 
respectively. The mean velocity components v$ and can be found by matching (6.7) 
and (6.8) with the outer solutions (3.2) and (3.3) as 

= = = -S' = 0. (6-24) 


where 6^ and are constants. The transverse velocity components vff and are non- 
zero because of the viscous spreading of the mean flow as mentioned in §3. The nonlinearly- 
generated Z-independent components of the spanwise velocity , W / ^oy[,*,m,n]» are not zero 
inside the critical layer although they vanish on its edge. 

If we substitute the expansions (6.10), (6.14) - (6.16) for u^, (6.11) for w (1) and similar 
expansions for v^ k \ and along with the mean-flow components (6.12), (6.21), (6.22) 
and (6.24) into the nonlinear equations in Appendix A and if we use the relation 

min( 6 

£ F -* Ci ' (6 - 25) 
j|o /|6 j|a+6 ^=max(— 6J,— aJ±j) 

where the summation notation defined by (6.19) is used and a and b are positive integers, 
then we can obtain a system of partial differential critical-layer equations, which will be 
shown in the following section, after considerable manipulation. The equations were given 
in Goldstein & Lee (1992) and Lee (1997a) for the single- resonant-triad interaction. The 
system of critical-layer equations can be solved analytically (Goldstein & Choi 1989; Gold- 
stein & Lee 1992, 1993; Wu 1992, 1995; Wu et al. 1993) or numerically (Lee 1997a) with 
the transverse boundary conditions at fj = ±oo that can be found from the matching with 
the outer solutions (5.30) to (5.33). 

The outer solutions (5.30) to (5.33) are obtained for the case when the oblique-mode 
amplitudes with positive and negative propagation angles are the same. We can show from 
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the critical-layer equations and their boundary conditions that 

l J a,b\3[,t,m,n] ^ ' a,- 


-fry[,f,m,n] 


(6.26) 


for all daggered components, i.e. ()*, that appear in the expansions (6.14) to (6.16) and 




s.1 (2)* 


(6.27) 


where the asterisk denotes the complex conjugate. 


The solutions of the critical-layer equations have to match onto the discontinuous 
0{6zda r ) = 0(cr 4r+1 ) and 0(e 2 d^ r ) = 0(cr 5r+1 ) terms in (5.30) and also 0(6 3c icr T+2 ) and 
0((2d<7 r+2 ) terms in (5.31) for the oblique and two-dimensional modes, respectively. It 
follows from (5.30), (5.31), (6.6) and (6.7) along with (6.1) and (6.5) that 

f v^dfj = -|idt f u^dfj = Re2/i c 7y c - <^j) A } E CJ + • • •> (6.28) 

7—00 J — OO J= — J 

r x )jgdil = -id r u^dfj = £ R efi c aY c - 4>j) (6.29) 

J — OO J — OO J= — J 

where E 3 and E CJ are defined in (5.34). The phase jumps in (6.28) and (6.29) can be 
eliminated by using (5.27) and (5.28) to show that 

Ijgdri* £ ( cos9+ ^«) f (^t) a ‘. e ‘” 

(6.30) 

£ 4 % .£*%£[*,- + 1 (f ) 1/2 } *\ («» 

From the expansion for t>( 3 ) that is similar to (6.15) and the expansion (6.16) for 
along with (6.26) we can show that 

j r 

—(3) Vr.of ^( 3 L) ^(3a) + + y( 3 «) + V (v {3d K , 

V y 1 = 2^ Re2£ C j V lly + V 1A . 0j + V^ U} + yi,ly,0+ *T,ly,0 + 2^ l/l,ly/y-* 

, 0 7 i * 




+ V i,v,],e,j + + 2^ v i,iyjtj +•••« 

e=-J~] J 
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u (4 > = Re2 ^ 
3=-J 


J +J 

i fr( 4a ) i r/( 4b ) 4. r/( 4<i ) 4. //I 4 *) 4. V'' ( //( 4c ) 

^2,0y + t '2,0;2j,j + ( '2,0y,0 + U 2,0; } ,0 + U 2,0;2j,] + 2-r V^ 2 - 0 ^ 

J-j ^+2j 

, r>( 4e ) 4//W) 1 iy( 4 s) \ 4. V 77 < 4fc > + V If7 (4j) 

+ t/ 2,0;2j,O + ^ 2, O', 2j, 1,1— j + U 2,0;2],t,]-t) + 2^ U 2,0',2j,t,j + 2^ \ ty 2,0;2j,/,<-J 

(=-J-j t=-J+ 2] 

min( J+j f J+£) 

E 

m=max(- J+j,- J+£) 

min(J— j, J— £) 

E 

m=max( — J— j,— J— t) 


1 r 7 ( 4n ) 4. 

+ t '2,0;2j,O-/ + 


( Tj ( 4k ) 4- /7( 4 ') 4- 7 /( 4m ) ) 

\^2,0;2j,f,mj-tn ' u 2,0;2j,t,m,} • u 2,0',2j,(,m,m-j) 

(4°) | 


' r y - ^2,0;2j,t,m,j 




(6.33) 


where we have only shown the components of the base subharmonic and fundamental fre- 
quencies. Substituting the above equations into (6.30) and (6.31) and by equating the 
coefficients of E C] and E } for the oblique and plane waves, respectively, we obtain the jump 
equations that will be given in the following section. They arise from the requirement that 
the velocity jumps across the critical layer which are calculated from the internal critical- 
layer solutions axe equal to those calculated from the external solutions. The jump equations 
must be solved with the system of critical-layer equations to determine the unknown insta- 
bility wave amplitudes A 0 and B : . 


7. The system of critical-layer equations and the jump equations 

In this section we will present the system of partial differential critical-layer equations 
along with the transverse boundary conditions and the jump equations. These final equa- 
tions are normalized in such a way that their nonlinear growth parts are free from the mean 
flow dependent parameters (Lee 1997a). 

It is convenient to introduce the following normalized variables (Goldstein & Lee 1992; 
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Lee 1997a) 


x = k ^r^axi - x 0 ) , V - v/(kc) ~ 2 ^oj X = X - 


2X 


k 3 r w ac 3 ’ 


{x,x> = - — {x,x}, 


and 


i( 2 ) - 


2irY r a 


KT U 


*Y c a 


°1 M ~ £-3 


c« (2) 

~T a \M - TZ 3=1*0 


(7.1) 

(7.2) 

(7.3) 


where the positive real constant k is a normalization parameter which can be chosen arbi- 
trarily, x 0 , rj 0 and Xo are coordinate origin shifts to be chosen subsequently, x\, X, X and fj 
are defined by (2.3), (3.5), (2.8) and (6.1), t w is given by (3.15) or (3.16) and x is given by 
(5.23). The frequency-detuning factor x, introduced by (2.6), can be given as an input for 
a problem. The mean flow coefficient was defined by (6.21) and (6.23). As in before 
the subscript 0 has been omitted for the quantities of the Oth resonant-triad ( j = 0). 

Following Lee (1997a) the dependent variables are normalized as, for —J ^ j ^ J, 

A > = k 3 c 2>lyAj ’ - k 2 cs\n6 Vv ' 3Ql ’ (7,4) 

where Aj and B : are related to Aj and B 3 in (3.1) to (3.4) by (5.26), is the mode shape 
of the leading-order solution given by (6.10) and (6.11), 9 defined in (5.29) is 7r/3 at the 
leading order and we have put 


M = 8tt Yj? 


(7.5) 


V a .j(x) = exp [|ia {X 0 + r? 0 x ) - ij\x\ . (7.6) 

The other normalized variables {U, V, W, , n ] {U,V,W,P }^ ;j [ i<tTOin ], which are 

coefficients of the Fourier series given in (6.14) to (6.16) and similar ones for the transverse 
and spanwise velocities and pressure, are defined in Appendix B for k = 2, 3 and 4. 
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If we put 


„ _ d . - x d 2 
Cn ~dx +mT] x dv 2 ' 


(7.7) 


the leading-order equation of the system of critical-layer equations becomes 

Cl QW = A v (7.8) 

which can be solved with the transverse boundary condition 

QW - 0 as Tj — * ±oo. (7.9) 

The second-order equations are 

= -1 £ U & = -K L V& = 0, (T.10) 

— 1 

4><L = f » CqUqIj^ — 2lW^ 2 y £ + (^0^0, 2y/)^ ’ *0,2yA = 2 ^0,2y,<» 

(7.11) 


"l,ly — M.ly ’ 

*8? = (*? + H) + *!?£!(*)’ 


= -<5-( sill2(? ) [q^ + Kjx- 

^ 0 )Q5 1) ]+#ij(sin2<?)xgf) + A„ 

(7.12) 

v!i L J = B„ V <%> = wgg = 0 

(7.13) 


f£iW^ 2 ) ^ yl 2 l — 2U^ 

^2 **, 0 y,< y ^ » V 2 , 0 Ul Aj “ 

(7.14) 

= iE-V'C™. 

it( 2 ) _ _w( 2 ) y( 2 ) _ n 

^" 2 , 2 y ~ " 2 , 2 y’ V 2 , 2 y ~ u > 

(7.15) 




and 


(7.16) 


p( 2 ^) _ p( 2 ^) __ n p( 2 ) = n 

M.lOT ‘*2,0 \jn U > r aMj[r ’ 

where the summation notation is defined in (6.20), is given in Appendix C and 

v[ 2 ^(x) could be found, although the present analysis does not require that, by matching 
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with the outer solution (3.2), ( 5 . 2 ), (5.31), (6.1), (6.5) and (6.7). The subscript 77 denotes 
the differentiation with respect to rj, i.e. d/dr ?, as usual, and the subscript {a, 6 } in (7.16) 
denotes either {0,0}, {0,2}, {2,0} or {2,2}. We have also put 

= \Aj§^ + idr(sin 2 9)Q^ for d = 0,±1,±2,..., (7.17) 


and 

r = c!{t w Y c ), 


(7.18) 


where r becomes unity at the leading-order if we use (3.20). Some of the simple equations 
have been solved with appropriate transverse boundary conditions and we have also used, 
in (7.11) and (7.15), the solutions of (7.16), 


p(2) 

r o,2;j,t 



(7.19) 


From (3.4), (3.11), (5.2), (5.33) and (6.9), we can find but it does not play any 

major role in computing the amplitudes A : and B r 
The transverse boundary conditions are 


w( 2 ) 77 ( 2 ) r/( 2 ) 0 as r) 


±00. 


The solutions of (7.11) given in Part 2 show that = 00 ) - = 

v o%A' n = °°) ” ~ _0 °) are non-zero, therefore, 


77 ( 2 ) _» 77( 2 ) ±0 ° 

U 0,2;j,t u 0,2;j,e ’ 


^( 2 ) 

*0,2y,< 


[>( 2 )±<» 

y 0 , 2 y,/ 


as 


± 00 , 


(7.20) 
— 00 ) and 

(7.21) 


where and V^y ,? 0 are functions of x. 

The velocity jump across the critical layer for the oblique mode will be obtained by 
solving the following third-order equations 


r w( 3L ) _ _ p( 3i ) _ r r y( 3L 

“ M.ly "t ”i,ly > 


(3i) - 2 : fi (2)j 

Vyjrm - Zla l M a v 


(7.22) 
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= - {fg"i,.o} + {l,2^}(co (7.23) 

(36) _ 


r , wl°°J — _ d( 36) , i_ 


E W U - + E ***<», 

/:[2] <:[-2] 


y(36) _ q 

* K Uu u ’ 


r w (3<0 - p( 3d ) 4 . r-< 2 l°> 1:1 - 1 ’ 1 ) r y( 3< *> - f/^M^.o) 1 

"'l.lyAm - M.laAm + v l,ly,r,m»jTj - v > 

P w( 3e ) -_p( 3e ) , ^(2|0,l:-l,-l -3) A y(3e) = _ r^(2|l,-l:2,0,-2)l 

A ' 1 "l,laAm ~ r l,l;),(,m +U (-j+/)|2,0;/,m » A ' 1 l l,ly/,m»j>j [°(-J+<)|2,0^,m J , 


(7.24) 

(7.25) 

(7.26) 

(7.27) 


f w(3/) _ _ p( 3/) _ ^(210,1:1,-1,3) p y( 3 /) — _ (V^ 2 ! 1 ’ 1 =2,0,2) 1 (7 28^ 

X 'l IV l,laAm” M.l \jJL,m L ’(-J+/)|2,0^,m> M.lyAmijij “ L (“J+^^.O^mJ ’ { < -*0) 


p w(3fl) p(3s) _ \p... ,iy( 2 )* /’.vW — i A 

*■'1 ” l,ly,< •* l,ly,< 1,f (l)b+< ,r 0,0^» ^l,ly,/»ji) — ,/1 J+< l/ 0,0;f»)T)) 


(7.29) 


r w( 3/l ) _ p( 3fc ) 4. A(2|0,l:-l,-l,l)* r |/(3fc) _ _ f^(2|l,-l:2,-2,0)* 

^ W l,ly/,m ~ •*l,ly,/,m **" U (j+/)|0,2;<,n» » A ' 1 K l,ly,*,m»j7j |^(j+/)|0,2;/,m J ^ ’ 


(7.30) 


with 


<,^rC*HrC[/,l + R S:; 3 W fM fc = 3i,3.,...,3A, (7.31) 


— 1/4*) 


7 ( 6 ) 


,w 


and 


?(*■) 


l,ly[,/,m] 7 j “ 0 ^° r * _ 3 ^’ 3<I ’ ' • • ’ 3 ^> (7.32) 

where ; , in (7.31) is equal to zero if k ^ 3L and 72.^] is given in Appendix C along 

with W[ in (7.22), We have put, in addition to (5.29), (6.20), (7.3), (7.17) and (7.18), 
that 

{^ZZT - * (*°s« + 

+*i(sin 2 9 ) {e! 1 ’, , ( 7 . 33 ) 
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where d\ to d$ are real constants. The pressure Pj j y j ( m j(x) for k = 3 L, 3a, . . 3/i could 
be found by matching with the outer solution, but it is unnecessary in this analysis. 

The above partial differential equations can be solved subject to the transverse bound- 


ary conditions 


l"i,ia + "l.ia v i,i \prn' v i,ialA* 


0 as rj — *• ±oo, 


(7.34) 


for k = 3a, 36, . . 3 h. 

The other third-order equations that must be solved in order to compute the velocity 
jump of the plane wave across the critical layer are 


V%J = B„ V,«“» = 2 + «©*). <g = °. (7.35) 


cMla = V,® A = 2V<%j + 

£3 = - {^,0} + (cos 2 «)vS' {£,-2$,} 

UVj^, = -2i > 

r w^ 3c ) — p( 3c ) 4. ^(2|0,1 : -1,1,-1) /'oV’( 3 c ) — f/7( 2 l 3 ' 1:2 > 4 - 2 ) 1 

^ 3 " 31y < m - T ’ ^3 v 3,l;],e,mnv ~ L (j— ^ ’ 

r U/W) _ p(3d) _ ^,(210,1:1,1,1) _ [^(213-1:2,4,-2)] 

L 3 W Z,V, J< e,m ~ _/ 3 fl y,f,m U (j-f)|2,0;r,m> ^3 ~ J ^ ’ 


(7.36) 


(7.37) 


(7.38) 


(7.39) 


(7.40) 


U 3%,([,m] - 3 (^3,ia,/[,m]»j “ for k ~ 3a, 36, 3c, 3d, (7.41) 


^ ( S = ^oU = 0, P 3 ( 5U m] , = 0 for k = 3a, 36, 3c, 3d, 


(7.42) 
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where we have used (5.29), (7.17) and (7.33), and in ( 7 - 36 ) are S iven 

in Appendix C and t4foy($) can found fr° m ( 3 - 2 )> ( 5 - 2 )> ( 5 -31)> (6-1), (6-5) an( * (6.7) 
although we do not need to obtain it. 

The transverse boundary conditions are 


w: 


(3) 


2,0-,j,O 


TT (3) wW 

u 2fl-,j,tv’ Kr 3,ly4,m]> K 3,l m]n»? 


0 as 77 


± 00 , 


(7.43) 


for k = 3a, 36, 3c and 3d. 

Finally, the velocity jump across the critical layer for the two-dimensional instability 
wave can be found by solving the following fourth-order equations: 


CjVlS, = 4i«S,i»„ 


(7.44) 


C.P&U = c 2 u!$l Jn = = i Offl,' < 7 - 45 > 

= £- 1 ^ 2 - -^(tan 2 0) (^O.Oy-i + 7 0,0;(-j+<)) ’ 

(7.46) 

A <U<,4 - i(sec 2 »)* • <™> 

- (-c 2 «) £>C» = (**’«> < 7 - 48 > 

Z’ Tj(4c/4j) _ / 2 m [/^(3o/3e|l,0:l,3,-l) -O|(3a/36|l,0:l,l,-l)l 17 49') 

n rr(4//4m) _/• 2 n'i f/~<(3e/3/|l, 0:1,3, -X) A(3c/3d|l, 0:1, 1,-1)) 17 501 

£ 2^2,0y,/,m,nn “ l SeC [ Cr (j-<)|l,l/,m,n ^(-j+^la.l^m.n J„ ’ l'- 00 - 1 

f, ( rr(4n) rr(4o) \ _ (' c „.2 /i\ 9 f z~<( 33|1,0:1,3,— 1) ^-,(3/i|l, 0:1,3,— 1) \ 

^2 {02,Oy,«,m»»> ^2,0 y/.m.njj / ~ ( SeC \ G (j-*)|l,l;<,m ’ C7 (j-/)|l,l^,m,nJ > V ‘ * &1 ) 

and 


V&Amj*, = f ° r 1 = 4i - 40 4<> ' 


- ottW 


,(k) 


(7.52) 
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where in ( 7 - 45 ) and in ( 7 - 52 ) are S iven in Appendix C and D^o^.m.n] = 0 if 


k ^ 4a. We have put 


rr(k\dud 2 ) 

n a,b;t,m 


^e)(i,u[ k l Am + i 2 w[ 


( k ) 


- V' 


(k) 


a ,b;£,m v a,b;t,m 


A.) 

1 dr}) ’ 


(7.53) 


and we have also used (5.29), (7.3), (7.18) and (7.33). The equations for and 

U-zo-j t mr, are obtained from (7.49) by choosing the superscripts {3a, 3a} and {3c, 36}, re- 
spectively, and similarly the superscripts {3e,3c} are used for U^. } im nT) and {3/, 3d} are 
for in (7.50). 

The solutions of the above fourth-order equations must satisfy the following transverse 


boundary conditions: 


Tj(k) 

U 2,0;j[,t,m,n]ri 


0 as tj — *■ ± 00 , 


(7.54) 


for k = 4 L, 4a, . . ., 4o. 

The above system of partial differential equations must be solved with the jump equa- 
tions in order to obtain the streamwise evolution of A : and B y The jump equation for the 
oblique mode becomes from (6.30) and (6.32) 


(cosd+ jij) | Ajt+ jOx - 1«Mj] - (2a r kT w X)' n A, 

= - •) /_“ *» [V, S3» + vSim + v uL + V&v, + F S.o,„ 

, V" ( y(3d) , v ( 3e) , ^(3/) , y(3h) \ 


and the plane-wave jump equation becomes from (6.31) and (6.33) 


(7.55) 


B„ + iOx - -fc) B, - |f B, = -±- f°° dr, [| V%>„ + U^,,„ + CU 

71 T J —oo 


J+j 


rr (4d) TT {4i) ST fr;( 4c ) j. M 4e ) 4 . r/( 4 /) + TJ^) \ 

+ U 2,0;j,0 V + U 2fi-,2},m + 2^ + V 2,0;2],t,n + U 2,0;2j,t,(t-j)v + U 2 , 0 ; 2 ^,(j-^)r j ) 

e=-J+: 
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t=-J-3 

+ E 

m = max ( — J+j,— J +£) 


J-3 J+l3 

-u V jj^ h) _l v / + t/ (4n) 

+ 2 ^ U 2fi;23,(,3n + 2_y \ V 2,0;23,t,(t-3)v + U 2,0;2 3,(,(3-t)n 


e=-J+2 } 


( tt^ 4 ^ + u {4l) 4- r/( 4m ) ^ 

min ( J — j , J — £) 


+ E 


M io ) 

2 , 0;2 3,t,m,3V 


, (7.56) 


where (5.23), (5.29), (7.1), (7.2), (7.4), (B3) and (B6) have been used along with 


t = a/(Twc), (7.57) 

and f (also in (C 11)) becomes unity, at the leading order, if we use (2.5), (5.8) and (5.9). The 
terms involving ijx on the left hand sides of (7.55) and (7.56) account for the wavenumber 
difference between the jth and Oth resonant-triads (see (5.21)) which is caused by the 
frequency detuning introduced in (2.6). In this analysis, it is convenient to choose the 
transverse coordinate origin shift T)o introduced by (7.1) to be zero 


Vo = 0. 


(7.58) 


The upstream flow can start out as a system of resonant-triads of linear instability waves 
(Goldstein & Lee 1992; Lee 1997a), as a system of resonant-triads of linear and nonlinear 
instability waves (Wundrow et al. 1994; Goldstein 1994; Lee 1997a) or as a system of pairs 
of nonlinear oblique modes (Wu et al. 1997). The upstream boundary condition in the 
linear instability-wave case becomes 


Aj -+ a, exp [(k 0 6 - \jx) x > Bj -*• b 3 exp [(* 2 d - i?x)*] as x -°°> ( 7 - 59 ) 


where 


** Ur/ 


1/2 / ^4 \ K ob/^2d”l/2 p(0) 

( K * \ / : a - \ 7 /. - - x 

exp JKXX 0 J , bj = exp (ijk X x 0 ) , 


MW 


d(°)| 


(7.60) 
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X and M are given by (7.2) and (7.5) and we have used (5.23). The A\ 0 ^ and i?j°' are the 
complex initial values of the oblique and plane wave amplitudes A 3 and B 3 in (3.1) to (3.4), 

Aj -*• A< 0) exp , B 3 -* 5j 0) exp (jkT w aKuXi) as x x -»• -oo. (7.61) 

The origin shifts i 0 and X 0 in (7.1) are chosen to satisfy 

(M/k 4 )B ( Q ) e iXo+kK2dXo = 1. (7.62) 


The upstream matching conditions for a system of resonant-triads in the later stage of 
the critical-layer evolution can be obtained from the solutions of the preceding linear or 
nonlinear stages as in Wundrow et al. (1994), Goldstein (1994), Wu et al. (1997) and Lee 


(1997a). 

The linear growth rates of the oblique and plane waves, k 0 (, and «2 d, can be obtained 
from (7.55) and (7.56), 


.1/2 


(7.63) 

(7.64) 


** = -ds; 

The last terms in both equations are due to the viscous Stokes-layer effect and this effect 
enhances the linear growth rates of the instability waves. Their magnitudes are of 0(cr r/2 ), 
but they will become 0(1) in the viscous limit where A = 0(a~ r ). 

It is easy to solve (7.22) and (7.44) to show that, as in Part 2, 

/” = 2br«S W [" iT[V (7-65) 

J — oo J —oo 


where is given by (7.3). As mentioned earlier, see (6.21), (6.23) and (7.3), (that 
is composed of two parts) is proportional to the double derivative of the base mean flow at 
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the critical level. The first term on the right hand side of (3.19) is due to the mean pressure 
gradient and the adverse pressure gradient (positive value of /2) increases the linear growth 
rate, while the second term, which is always negative, is due to the Blasius profile. In the 
later non-equilibrium critical-layer stage the exponent r of the growth-rate parameter a T is 
smaller than 3 and also /i < 1 as shown by Wundrow et al. (1994), Goldstein (1994), Wu 
et al. (1997) and Lee (1997a). Thus, from (3.19), the linear growth rates given by (7.63) 
and (7.64) become negligibly small compared to the nonlinear growth rates. 

8. Concluding remarks 

The non-equilibrium critical-layer analysis of a system of frequency-detuned resonant- 
triads is presented. It is shown that resonant-triads can interact nonlinearly between them- 
selves in the common critical layer if the scaled Strouhal numbers (of their fundamental 
plane waves) are different by the magnitude of the growth-rate parameter, 0(a T ). The flow 
outside the critical layer is still governed by the linear dynamics. The generalized scaling 
of Lee (1997a) has been used for the long- wavelength small-growth- rate instability modes 
in boundary layers with and without mean pressure gradient. The wavenumber parameter 
a characterizes the small wavenumber and the growth-rate parameter a T characterizes the 
ratio of the small growth rate to the wavenumber of the instability wave. Although the anal- 
ysis was carried out with the non-equilibrium critical-layer scaling, the quasi-equilibrium 
amplitude equations can be obtained in Part 2 by taking the viscous limit of the finite- 
viscosity amplitude equations. It was shown by Goldstein (1994), Wu et al. (1997) and Lee 
(1997a) that the later downstream stage of the quasi-equilibrium critical- layer evolution is 
governed by the non-equilibrium dynamics. In fact, the final critical-layer stage will be 
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governed by the inviscid version of the system of equations presented in §7 without the 
linear growth terms. Since the value of r becomes nearly equal to zero there, the viscous 
effect and the linear growth rate will become negligibly small (Wundrow et al. 1994). 

The effect of the viscous Stokes layer has been included so that the linear growth rates 
are valid for both in the finite-viscosity and viscous-limit cases. The upstream instability 
waves can be composed of either linear instability waves or nonlinear waves of the preceding 
critical-layer stage. 

The nonlinear effect in the critical layer produces a term which is periodic in the span- 
wise direction with the spanwise wavenumber equal to twice of that of the oblique modes. 
In contrast to the previous results (Goldstein & Choi 1989; Goldstein & Lee 1992, 1993; Wu 
1992, 1993; Leib & Lee 1995), this term becomes unsteady in the present multi-frequency- 
interaction case. The magnitude of this nonlinearly generated low-frequency mode is as 
large as that of the primary oblique modes outside the critical layer. The complete so- 
lutions of the span wise-periodic components, as details will be given in Part 2, can be 
obtained from the multi-layer analysis as in Wu (1993). When r = 3 and A defined by 
(2.8) is 0(1), we need to consider the solutions in the potential region where y = y/cr , 
the main boundary layer where y = 0(1) and the viscous wall layer where y = <t 3 Y in 
addition to the inviscid wall layer of 0(cr) and the critical layer of 0(a 4 ). The viscous wall 
layer, which is thicker than the viscous Stokes layer in §4, is required in order to satisfy the 
wall boundary conditions. (The inviscid- wall-layer solutions of the stream wise and spanwise 
components become singular at the wall.) The three-dimensional boundary-layer equations 
in the viscous wall layer must be solved with the no-slip boundary conditions at the wall and 
appropriate boundary conditions on the upper edge of the layer in order to match with the 


45 



solutions in the inviscid wall layer. The transverse velocity on the upper edge of the viscous 
wall layer is determined internally by the boundary-layer equations themselves and it can 
not be given as a boundary condition as in Wu (1993). It is interesting to find out that the 
velocity components in the inviscid wall layer, at the leading order, become equal to zero 
where Y > Y c but non-zero where Y < Y c (Lee 1997c). Thus, the nonlinear interaction 
between oblique modes produces, at the leading order, the spanwise-periodic velocities only 
below the critical layer. The non-zero velocities at Y = F c _ will be determined from the 
matching with the critical-layer jumps. 

Since the velocity jumps of these low-frequency modes across the critical layer are 
determined from the solutions of the amplitude equations, which become singular at a finite 
downstream position, the magnitude of the spanwise-periodic velocities also becomes very 
large near the singular point. When r becomes 9/10 the streamwise velocity becomes as big 
as the base mean flow, therefore, the nonlinearly generated spanwise-periodic component 
and the base mean flow start to interact nonlinearly in the viscous wall layer whose thickness 
is of 0(<r 37 / 10 ) (Part 2; Lee 1997c). This nonlinear viscous wall layer will be separated into 
two layers (inviscid nonlinear layer and the viscous nonlinear layer) in the later downstream 
region. 

When the difference between the scaled Strouhal numbers, s } given by (2.5), of the 
resonant-triads is larger than the magnitude of the growth-rate parameter <x r , the instability 
waves of an individual resonant-triad initially grow independently of the instability modes 
of the other resonant-triads. For example, when r = 3 in the adverse-pressure-gradient 
boundary layer (Goldstein & Lee 1992), the plane and oblique waves of the resonant-triad I 
(and those of the resonant-triad II) go through the usual linear stage, parametric-resonance 
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stage and the fully coupled stage, independently of the other resonant-triad, as shown in 
figure 2, provided the (fundamental) frequency of the resonant-triad II is larger than that 
of the resonant- triad I by 0(a p+2 ), where p is a positive number that is smaller than 3. 
Since the distance between the critical levels, j/ c // — y c i , which is of 0(a p+1 ), is larger than 
the critical-layer thickness of 0(cr 4 ), the growth of each resonant-triad is not affected by 
the presence of the other resonant-triad. The streamwise evolution of the oblique and plane 
waves are still determined by the integro-differential amplitude equations (5.50) and (5.51) 
of Goldstein & Lee (1992), or by the 0(1)- viscosity equations given by Wu (1995) and Lee 
(1997a). 

The non-equilibrium critical layer becomes thicker and r becomes smaller as the insta- 
bility waves propagate downstream (Goldstein & Lee 1992; Wundrow et al. 1994; Goldstein 
1994; Wu et al. 1997; Lee 1997a) mainly due to the self-interaction between oblique modes. 
The instability modes of different resonant-triads start to interact nonlinearly within their 
overlapped critical layers when r becomes equal to p. The oblique and plane waves of both 
resonant-triads I and II in figure 2 become fully coupled in the frequency- detuned nonlinear- 
interaction stage where their growth rates are now 0(a p+1 ), or r = p(< 3), which is the 
initial difference between the respective critical levels. The amplitudes are then determined 
by the inviscid version of the system of equations presented in §7. When r approaches to 
zero the critical layer becomes almost as thick as the inviscid wall layer whose thickness is 
of 0(a) and, thus, a wide range of instability waves with the scaled frequencies differing 
by alm ost. 0(1) can now nonlinearly interact between themselves. This non-equilibrium 
critical-layer stage is still governed by the parabolic type equations of this paper, but the 
flow in the next stage will be governed by fully elliptic triple-deck equations (Goldstein & 
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Lee 1992, 1993; Wu et al. 1997). The turbulent flow field may be determined by solving 
these elliptic equations with the upstream matching boundary conditions that can be ob- 
tained from the solutions of the equations presented in §7 (with appropriate downstream 
boundary conditions). 

The system of critical-layer equations are presented in §7. The variables are normalized 

* 2 ) 

(Lee 1997a) in such a way that the mean-flow-dependent parameters, i.e. and t w , only 
appear in the linear growth rate terms in (7.63) and (7.64) or in (7.22), (7.44), (7.55) and 
(7.56). We can show that a^, and b$ in (Cl), (C2), (C4), (C5) and (C7) 
do not play any major role in computing the velocity jump across the critical layer and 
they will not appear in the final amplitude equations (Goldstein & Lee 1992; Lee 1997a). 
The no nlin ear part of the system of critical-layer equations is free from any mean-flow- 
dependent parameters except A. Therefore, the results of this frequency-detuned analysis 
can be universally applied to other shear flows, with minor modifications for the specific 
linear growth rates, which have been individually studied for the single resonant-triad or 
single- frequency oblique modes, for example, oscillatory Stokes flows by Wu (1992, 1995) and 
Wu et al. (1993), free shear layers by Goldstein k Choi (1989) and Mallier k Maslowe (1994) 
and supersonic boundary layers by Leib fe Lee (1995). The magnitude of the wavenumber 
of the instability wave in these analyses is of order one and thus the near-neutral small- 
growth- rate approximation that was originally introduced by Goldstein k Leib (1988) and 
Goldstein k Hultgren (1988) has to be used instead of the long-wavelength small-growth- 
rate approximation used in this analysis. The nonlinear part of the system of critical-layer 
equations is also valid for any values of 0 (for the oblique-mode interactions) since the 
resonance relation (5.15) or (5.16) has not been used in it. 
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The frequencies of the fundamental plane waves of the frequency-detuned resonant- 
triads can be freely specified (or can be given by experiments). The small parameter cr, 
which accounts for the small wavenumber of the instability waves in boundary-layer flows, 
can be chosen from the adverse-pressure-gradient boundary-layer scaling (Goldstein et al. 
1987; Goldstein & Lee 1992; Wundrow et al. 1994), the upper-branch scaling (Mankbadi 
et al. 1993; Goldstein 1994; Wu et al. 1997), or the favorable-pressure-gradient boundary- 
layer scaling (Wu 1993). The scaled Strouhal number s 3 and the frequency-detuning factor 
X are then computed from (2.5) and (2.6). The wavenumber of the plane wave a 3 that 
possess the series expansion (5.8) is fully determined by the dispersion relations (5.9) to 
(5.13) and the phase speed c : can be obtained from (2.5). The wavenumbers ji } and 7j 
of the oblique modes are computed from (5.15), (5.16) and (5.18). In order to obtain the 
streamwise evolution of the (complex) oblique and plane wave amplitudes A : and B } of 
the full system of resonant-triads, the system of partial differential critical-layer equations 
in §7 must be solved with the jump equations (7.55) and (7.56) along with the upstream 
conditions (7.59). However, we only need to solve (7.8) - (7.16), (7.20), (7.22) - (7.32) and 
(7.34) with (7.55) and (7.59) for the nonlinear interaction between the frequency-detuned 
pairs of oblique modes. 

The system of partial differential critical-layer equations can be solved both analytically 
(Goldstein & Choi 1989; Goldstein & Lee 1992, 1993; Wu 1992, 1995; Wu et al. 1993) and 
numerically (Lee 1997a). It may be more convenient to solve them numerically for the 
finite-viscosity resonant-triad case. But it is relatively simple to solve them analytically 
when only the pairs of oblique modes are considered. The resulting analytical amplitude 
equations, as will be given in Part 2, are quite compact and can be easily solved numerically. 
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Comparison with the analyses of Goldstein & Lee (1992) and Lee (1997a) shows that the 
previously considered single- resonant-triad is a special case of the present multi- resonant- 
triads. The perfectly tuned case of Goldstein & Lee (1992) and Lee (1997a), when k, in (8.3) 
of Lee (1997a) is equal to zero, corresponds to the interaction of a single Oth resonant- triad, 
i.e. between the Oth plane wave and the pair of Oth oblique modes. 

The system of partial differential critical-layer equations in §7 will be solved analytically 
and the amplitude equations (without the back-reaction term in the plane-wave amplitude 
equation) and their viscous limit will be presented in Part 2 (Lee 1998b) along with the 
numerical solutions. The low-frequency spanwise-periodic components that are induced by 
the nonlinear effects in the critical layer will also be considered in Part 2. 

This work was supported by the Acoustics Branch at NASA Lewis Research Center, 
contract number NAS3-98022. 

Appendix A. Governing equations for t>^ f and p W 


au^ + = 0, (Al) 

au l J? + 4*- 1} + 4 fc) + = 0 for k = 2,3 and 4, (A 2) 

J 

Lu^ = ^ Re2ir u ,7c(sin 2 0)A J E CJ , (A3) 

j =-•/ 

J 

TkjI 1 ) = ^2 ReT w ac(smO)A J E S] , (A 4) 

j=-J 

Lu? ] = 4V Y Ke2iT w Y c iA}E C3 + t w - fju + ^xatpW) 

j =-J 

- -1- , (A 5) 
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W 2) = [V3jxc{cos 2 0)ReAjE SJ + ^ 1) Re(2ir c 4£' CJ )] - - T w fjw^ 

—au^w^ — (A 6) 


j =-J 


Lu { ? ] = u*£} Re2ir w y c 7A J £ CJ + r w (fiw^ - fju ^} ,) + ^x« t w w^ - (w (1 ) u (j)) 


j=-J 


- [ti (1) U^ + 4 1) ^ 2) + 5 ( u(1)u(2) ) x +/ 3 ( ti,(1)u Z ) + ti,(2)u Z ) )]-’ ( A7 ) 

Lw {3) = J2 Re2ir w F c 7i JJ B CJ - h ( z ~ ^X« (4? + “ Ww™ 

J=~J 

— a + u^Wx^ — P ^ , (A 8) 


+^x<* 


j 

r 

3=~J 

T w w{ 


Lu£ ] = 4? S Re2ir w y c 7 J 4 J £ , CJ + t w (fiw^ - ifttJS) 
r 

ug) _ ^ ti ,( 1 ) tl g) -|- u,( 2 ) u g)) j _ [(ud)^)^ + «f)«P) + 4 2) «d) 

+Q ^ud) U ( 3 ) _|_ I u (2) u ( 2 )^ + ^ _|_ ^(2)^) + «?d) u d)^ ? (A 9) 


and 


p d* = 0 for k = 2 and 3, 


(A 10) 


where we have put 


r / d _ d \ __ d . d* /A tl v 

I = c fe + xa 0 tJ + rwa Vx- (A } 

E CJ and E 3J are defined by (5.34), x and X a 16 related by (5.23) and 6 is defined in (5.29). 
We have also used (3.8) and (3.20). 

Appendix B. Variables in the normalized critical-layer equations in §7 


J2£) _ 

01 , iy “ 


_ (Y c M) l/2 cos 0 V _ (2L) J 2 L) = c_M 


k z T w c 






(Bl) 
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{9o?0;j» ^Oy,*’ 9^1 y} ~ -3 tan 2 Q { V ^Aofi-p ^oAoJijp ^ A 2 fi-.jp ^y^Jy} > ( B 2 ) 

{^^Am]} = f ° r * = 3I ’ 36 — 3/l ’ 


( 3 L)_ M T) ^ _( 3 ) _ M C 0 f_e 43 ) rR4 x 

92 t 0;j “ TJ— ^2^92,0;^ ?2,0;j,* ~ £4 r ®2tf9 2l 0tf/» ' B4 ' 

rv / ^ ** * u/ 

f (3a) (36,3c, 3dn cM 3/2 COS 0 f -(3a) -(3b, 3c, 3d) \ fX) ^ 

Rli.O «S.ltf4.m]) = { P 3y + /9 3 ,ly/> ^^m]) ’ (B 5) 

f (4L) (4M (4c) (4d) (fc) 1 2zaM f 44L) -n -(46) 

{92,Oy»92,Oy,/’92,Oy,<,m’92,Oy,<»92,Oy,/l,TO,n]J — «5 r 3g \^2;2j02,Oy > ^2;2j-/92,0y^» 

^2y+m^0y,/,m> 7>2;2j+<?2,0y,*’ ^2y92,Oy,<[,m,n]} ^° r * = 4a, 4e, . . . , 4o, (B 6) 

where M and 2? a;j (i) are given in (7.5) and (7.6) and 9 is defined by (5.29). The normal- 
ization pair {q,q} of for k = 2L,2 ,.. ., 4o denotes either {U, U}, 

{V,2iF/(Kac)} 5 {W, W tan#} or {P,iP(tan 2 9)/(kt w c)}. 

Appendix C. w[™ ] in (7.12); w[% ] in (7.22); ttgg in (7.31); 2^, ftgL* 
and yl % } ,t i« (7-36); K^ t in (7.45); 3^ in (7.52) 


°)VW + jg t ) QW, (C 1) 

W}*"* = 2i7r£ (v ( t' X) vM + jg t ) [WJSS + 39 t 9 sQ^} + i *tg, [(v { / hl) ) 2 [ (fjV^Q^)^ 

+ |iA (4 4- 1>?>) <?«,} - 2ajJ^ l ^)g«] + ir* [4>i 2) + ®3M (** - ift) 

+ a?i {flV& + iA (l + §X>i 12/5) ) &}] ( c 2) 

^SS = 2 xi^ (Vlf lfi) u{%} + JftWfig) 7 (C 3) 
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^-lo-Mjw^-ogissr 

- 9 , + §«!-< Kw!. - 2»-®i' :, ’ OI <?S ,) )] , (C4) 

kSL = 0 - 20<w.w^ + [eSiS'l , - 9.fiv { f v) Q?.L< 

+ }«. [{i<i> - 2 S , (4' !l,0) - «*) Pr”} fljli] n > (C 5) 

, = -ij, + 0 - M***8U . (C 6) 

*g& = -2xto [« (4» s ,x>^'°»<j<X, ( + 4” 2 ’ 2) eS^) - o - 2t)s, (4 i'fcwjjSU 
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= -2xi« 9l K* 2 ' 0 ^, + 0 - *)*«&,]„ • (C 8) 

where A, «<£ , r, and jrgg ,1 - 1) are given in (7.2), (7.3), (7.18), (7.33) and 

(7.53). We have also put 

-pO- a ’ b ) = ■§= + i (jx “ \ ar to ) + i ba \] , + jiaA^y, (C 9) 


t) = Tj + \rj 0 , x = f + K®, 


(CIO) 


£ = a/(r w c), g t = ^ix cot 0, p s = sin 2 6, g e = sec 2 0, g n = 2i tan 2 0, (C 11) 
where ( becomes unity, at the leading order, if we use (2.5), (5.8) and (5.9). The mean flow 


coefficients are normalized by 


=(i) _ £al h (2) s (2) h {2) \ = 27rYca J a {2) 2q 3~m 2b M l (C 12) 

M kr 2 ' *• 2A *’ 3A *’ ft 2 r 2 c 2 ( 2M ' k 2 r w ac' k 2 ac 2 ) 

where a$, a^, and 6$ are defined in (6.12), (6.21) and (6.24) and a$ is given in 

(6.13). 
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Fig. 1 Multi-layer structure for the non-equilibrium critical-layer analysis of multi-resonant- 
triads in boundary layers. I, main boundary layer; II, inviscid Tollmien wall layer; III, 
critical layer; IV, viscous Stokes layer. 
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Fig. 2 Nonlinear interaction between the frequency-detuned resonant-triads I and II in an 
adverse-pressure-gradient boundary layer where p < 3. A, linear stage; B, parametric- 
resonance stage; C, fully-coupled stage; D, frequency-detuned nonlinear-interaction 
stage; E, triple-deck stage. 
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